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Abstract. We study experimentally the nonequilibrium phase behaviour of a 
horizontal monolayer of macroscopic rods. The motion of the rods in two dimensions 
is driven by vibrations in the vertical direction. Aside from the control variables 
of packing fraction and aspect ratio that are typically explored in molecular liquid 
crystalline systems, due to the macroscopic size of the particles we are also able to 
investigate the effect of the precise shape of the particle on the steady states of this 
driven system. We find that the shape plays an important role in determining the 
nature of the orientational ordering at high packing fraction. Cylindrical particles show 
substantial tetratic correlations over a range of aspect ratios where spherocylinders 
have previously been shown 1 to undergo transitions between isotropic and nematic 
phases. Particles that are thinner at the ends (rolling pins or bails) show nematic 
ordering over the same range of aspect ratios, with a well-established nematic phase at 
large aspect ratio and a defect-ridden nematic state with large-scale swirling motion at 
small aspect ratios. Finally, long-grain, basmati rice, whose geometry is intermediate 
between the two shapes above, shows phases with strong indications of smectic order. 
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1. Introduction 

Our experiments study the orientational ordering of horizontal, quasi-two-dimensional 
assemblies of millimeter-sized, elongated particles driven into motion by external 
vibration in the vertical direction. Since the particles we study are rod-like and do 
not deform appreciably, the natural point of comparison is the phase behaviour of 
thermal equilibrium hard-rod systems in two dimensions. The experiments address 
some basic questions on how systems in driven, nonequilibrium steady states differ in 
their phase behaviour and dynamics from similar systems at thermal equilibrium: Can 
the nonequilibrium phase diagram be rationalized by a comparison to hard-rod Monte 
Carlo simulations IJ alone or must we invoke new physics? In thermal equilibrium, the 
ordering of hard-rod liquid crystals is determined solely by geometrical considerations of 
packing, which should apply equally in the molecular and macroscopic realms. On the 
other hand, thermal fluctuations are of no importance in the systems we report in this 
article, and all exploration of phase space is made possible only by means of external 
driving forces. 

The nonthermal nature of the system has several implications: well-established 
extremal principles involving free energy and entropy cannot be invoked, the sources of 
fluctuation and dissipation are distinct, and nonequilibrium fluctuations may either help 
or hinder certain kinds of order. 'Moving XY models' of flocking |21IS], for instance, have 
been shown to display true long-range order in two dimensions, unlike their equilibrium 
counterparts. Nonequihbrium steady states with strictly nematic order jlj are predicted 
to have anomalously large density fluctuations, and velocity autocorrelations that decay 
much slower than in corresponding thermal systems. There are also a number of 
papers El 13 E] on the dynamics of orientable systems that are motivated by the 
biologically important motor-microtubule system. These systems support a variety of 
nonequilibrium structures, in particular, rotating vortices, which have no equilibrium 
analogue. While the particles in these systems are polar in nature and therefore not 
directly comparable to our particles, they do show the possibilities of qualitatively new 
phases and defect behaviours in systems of self-driven orientable particles P ITU]. 

In the context of nonthermal systems like granular media, Edwards and coworkers 
have taken an entirely different theoretical approach: they propose a statistical 
mechanics for granular systems by constructing an ensemble that includes with equal 
probability all mechanically stable realizations of a system. The phase behaviour of a 
system of anisotropic grains has been considered JT] within this framework. The theory 
predicts, remarkably, that elongated particles do not undergo an isotropic-to-nematic 
transition as a function of density. 

There have also been some recent experiments and simulations on the packing 
of anisotropic granular particles under gravity. Monte Carlo ^21 simulations of ellipses 
falling in 2-dimensions found that they tend to orient with their long axes perpendicular 
to gravity. Quasi-2D experiments and simulations on cylinders showed similar 
results, with nematic ordering extending over two particle lengths. Villarruel et al. Jl] 
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studied the compaction of rods in 3-dimensions: rods were poured into a cylindrical 
container and compacted by 'tapping' i.e. applying isolated vertical jolts. While 
quantitative measures of orientational order were not presented, it was observed that 
repeated tapping caused the rods to orient vertically, first at the walls and then in 
bulk, into smectic-like layers. Blair et al. ^H] have studied vibrated systems of rods, 
contained in circular and annular geometries. At high packing fractions, the rods 
separated into regions where they are stacked horizontally, and regions in which rods 
came out of the horizontal plane and are only slightly tipped from the vertical; in 
the latter regions, vortices formed, rotating in a direction governed by the tilt of the 
rods. A phenomenological explanation for the formation and motions of the vortices 
was proposed in Ref. JHl- Although we do report on nonequilibrium dynamical effects 
such as these, our major preoccupation here is with static orientational ordering. Also, 
the role of gravity in our experiments is to be distinguished from those reported in Refs. 
[T^ ir^l IT^ : while gravity is very important in the kinematics of the particle motion it 
does not play the role of an ordering field as in ^21 1^ 1^ • The ordering in our system 
arises purely from interparticle interactions. Measurements of the isotropic-nematic 
transition in a quasi-2D system of vibrated rods have also been reported by The 
focus in that paper is on understanding the long-wavelength distortions of the director 
in terms of competition between wall alignment and bulk elasticity. Experiments very 
similar in intent to ours have been performed by Fraden and coworkers |18j . 

While there are no firmly-founded theoretical predictions to guide the effort of 
charting out a nonequilibrium phase diagram for elongated particles, simulations of 
hard rods in equilibrium provide a valuable point of reference. Frenkel and others 
02] have performed Monte Carlo simulations in two dimensions, on sphero cylinders 
of length L and diameter D. They find that nematic phases do not form for particle 
aspect ratio A = L/D < 7. Instead, they identify an isotropic phase at low number 
densities and a layered solid at high densities. For particles with A> 7 a. nematic phase 
intervenes at intermediate area fractions. In a grand canonical Monte Carlo study in two 
dimensions, Khandkar and Barma [201 confirm the quasi-long-range nature of nematic 
ordering in the "needle" (i.e., D = 0) limit, and demonstrate that imposed gradients in 
the chemical potential suppress director fluctuations. A recent Monte Carlo simulation 
PT] studies the ordering of rectangles with A = 2 and finds a transition from isotropic 
to tetratic order. In an earlier work, Veerman et al. j22] report the occurrence of a 
cubatic phase in a system of truncated spheres and more recently Borukhov et al. jSB] 
find cubatics in linker-filament models for the cytoskeleton. While we are able to make 
a detailed comparison to these equilibrium simulations, due to the macroscopic nature 
of the particles we are also able to explore a new axis in the study of hard-rod systems: 
we can make essentially arbitrary variations in the shape of the particles by selectively 
etching cylindrical rods. We find a surprising dependence on the precise shape of the 
particle - objects with similar aspect ratios but dissimilar shapes order differently. It is 
not possible to determine on the basis of our work alone whether this shape-sensitivity 
is a nonequilibrium effect, however, we hope to stimulate more detailed investigations 
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Table 1. Tabulated below are the details of the various particles we used in our 
experiments. 





Aspect Ratio Length (mm) 


Rice 

Short cylinders 
Short rolhng pins 

Long cyhnders 
Long rolhng pins 


4.0 ± 1.07 7.1 ± 0.44 
4.9 ± 0.42 3.7 ±0.11 
5.2 ± 0.59 4.6 ± 0.16 
12.6 ± 0.66 9.9 ± 0.12 
10.4 ± 1.45 9.2 ± 0.39 



of these effects in hard-particle simulations. 

Following this introductory section is an experimental section, |21 in which we 
describe the particles used in our study, the geometry of the cell, and the driving 
mechanism. The system is visualized by real-space images; the acquisition and analysis 
of the images is also outlined. A detailed description of our results is in section 
but we summarize these briefly here: Strictly cylindrical particles show substantial 
tetratic (that is, four-fold orientational) correlations at both small {A = 4.9) and large 
{A = 12.6) aspect ratio. The length scale of the correlations increases as the number 
density is increased; however, we do not see long-range order (LRO) over the range of 
area fractions and vibration amplitudes we explore. Particles that are thinner at their 
tips ( "rolling pins" or "bails" , see figure d) show nematic ordering at high area fractions 
in the same range of aspect ratios. For shorter particles {A = 5.2) we obtain a defect- 
ridden nematic state with large- scale swirling motion at the number densities we study 
in our experiment. For A = 10.4, we find at high densities, a phase with at least quasi- 
long-range nematic order (QLRO), that is, power-law decay of nematic correlations. We 
can induce transitions out of this nematic phase either by lowering number density or by 
increasing the intensity of the driving. Finally, long-grain, basmati rice, with an aspect 
ratio of ^ ~ 4, shows phases with strong indications of smectic order. We comment 
on this below in section EISl We also find intriguing dynamical phenomena, with clear 
signatures of nonequilibrium behaviour, as discussed at the end of section 01 

2. Experimental setup and procedure 

Our experiments involve vertically shaking a horizontal layer of rod-like particles. The 
rods were confined to a single layer between an aluminium substrate and a lid that fit 
over the particles leaving them only a little vertical room to move in. In this section 
we describe the particles, the geometry of the experimental cell, the driving mechanism 
and, finally, the methods of image acquisition and analysis. 
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Figure 1. The 5 kinds of particles we used in our experiments. White silhouettes 
are schematic outlines of the 2D projections of the particles. We call particles with 
a strictly rectangular profile 'cylinders' and the others 'rolling pins' or 'bails'. We 
define the aspect ratio as ^ = L/D where the diameter D is measured at the largest 
cylindrical cross-section. Far right: basmati rice grains. For details of size and 
polydispersity of particles see Table H 

2.1. Particle Dimensions 

We used five different kinds of particles in the experiment. Figure ^ shows images of 
a few of each. The particles labelled "cylinders" are cut from a length of enamelled 
copper wire. The process of cutting introduces only a small degree of polydispersity 
(less than 5%) in the length of the cylinders. The particles labelled "rolling pins" (or 
"bails" |21j) are produced by etching the tips of the cylinders in ferric chloride solution. 
Since there are variations in the local etch rate, and in the exact length of rod that is 
exposed to the etchant, the degree of polydispersity in the bails is somewhat greater. 
The rice was used as purchased. Table [T] summarises the dimensions of the particles 
used in the experiment. 

2.2. Geometry of the Cell 

The particles were confined in a shallow circular geometry with a 13 cm diameter, as 
shown in figure El Whenever necessary, it was possible to further confine the particles 
to a square region inscribed within the diameter of the cell. The base of the cell is made 
of aluminium, anodized black to enhance optical contrast with the particles. The top 
of the cylinder is made of perspex, and is placed so as to leave a vertical space H = 1 
mm above the base. Since the particles have a diameter D ~ 0.8 mm, this gives them 
about 25% of their thickness to explore in the vertical direction. Thus particles cannot 
pass over one another by moving out of plane. The vertical degree of motion allows 
the planar projection of the particles to overlap a little when they are being shaken; 
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Figure 2. Top and side view of base of circular sample cell. The height of the sample 
volmiie, H = 1.0 mm, and the diameter of the particles, D = 0.8 mm. 

the horizontal component of the distance of closest possible approach of two particles is 
slightly less than D. Thus, the vertical clearance H — D can be seen as a parameter that 
tunes the particle interaction potential from hard- core repulsion {H — D = 0) when the 
motion is strictly 2-dimensional, to slightly softened repulsion ioi H > D. The shaking 
amplitude Oo can then be seen as playing the role of temperature. 

Another important feature of the cell is the boundary condition presented to the 
particles at the edge of the system. When the boundary is smooth, particles tend to 
lie parallel to the boundary. This occurs both when the boundary is a square or a 
circle since in the latter case the radius of curvature is much larger than the particle 
length. The parallel alignment is a consequence of dynamical stability: when a rod 
aligned with the wall collides with it at one end and rotates out of alignment, then the 
next wall-collision will be with the other end of the rod, which will tend to restore the 
alignment. No restoring torque of this nature exists for perpendicular alignment with 
the wall. In some situations, we induce perpendicular alignment by gluing a layer of rods 
normal to the wall. The tips of particles perpendicular to the wall tend to get trapped 
in the interstices between the tips of the glued particles. We found this procedure to be 
effective only for the bails, and not for the cylinders. 

2.3. Excitation Mechanism 

The cells were filled with particles to a known area fraction and then mounted on a 
permanent magnet shaker (LDS 406). For most of the experiments described here. 
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the cell was shaken at a frequency / = 200Hz and an amplitude resulting in a 
fixed peak acceleration, as measured by an accelerometer (PCB Piezotronics 352B02), 
of r = (27r/)^ao = 6g where g is the acceleration due to gravity. The corresponding 
amplitude typically used was Uo = 0.037 mm = 0.0461? , smaller than the clearance of 
H — D = 0.2D between the top of the particle and the lid. However, the characteristic 
velocity of the plate (27r/ao) is sufficient to launch particles to a height of ~ 0.1mm . In 
one instance we varied the acceleration between 5g and 8g to observe the effect it had 
on the degree of ordering. 

2.4- Image acquisition and analysis 

Dynamical features of the steady states were imaged with a high-speed video camera, 
however, for the majority of the results reported in this article, our primary interest was 
the instantaneous configurations of the particles. These were imaged using a consumer 
digital camera (Canon PowerShot G5) placed above the cell. Images were taken at a 
resolution of 2592 x 1944 pixels, at time intervals of approximately 15 sees. This interval 
was long enough for significant rearrangements of the particles to have taken place, and 
hence for each image to be a statistically independent sampling of the steady state 
arrangement of the particles. In the rest of this article, we will discuss data extracted 
from images that were taken while the sample was being continuously agitated by the 
shaker. Images taken by halting the drive (as was done in 18J) showed qualitatively 
similar behaviour but with noticeably greater disorder in the particle arrangement. 

The images were analyzed within the ImageJ j^H] analysis software. The images 
were subjected to a bandpass filter to remove pixel noise and long-wavelength variations 
in illumination, and then thresholded before identifying the centre of mass and 
orientation of each particle. The error in assigning the coordinates was about 0.5D 
(it was possible to reduce this error to 0.2D for the short cylinders) and the error in 
assigning the angles was less than 2.5°. Figure El has bars indicating the centres of 
mass and orientation of the particles superimposed on the image. Typically 97% of the 
particles were successfully identified, the remaining particles being difficult to resolve 
as independent objects because of the high packing density, and the consequent lack 
of intensity gradient between particles. The other objects identified are mainly particle 
pairs; we are able to determine the position and orientation of some, but not all, by 
using information on the minor and major axes of the composite object. However, due 
to the relatively small number of such instances, they did not affect the results we report 
here. 

3. Results 

This section is organized in the following manner: the first half will be devoted to results 
obtained using the cylinders and the latter half to those with the rolling pins. 

Figure El shows the region of phase space explored and gives us an overview of 
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Figure 3. Particles identified by image analysis routine: each identified particle was 
assigned a centre-of-mass coordinate and an orientation, indicated by the location and 
orientation of the bars in the image. Even at these high area fractions, typically 97% 
of particles were identified. 




Figure 4. Real space images of steady-states observed: striped, smectic-like state 
with the basmati grains (left), state formed by aspect ratio A — 5.2 cylinders with 
strong four-fold correlations (centre) and nematic state formed by ^ = 10.4 rolling 
pins. 

the results. The cyhnders, irrespective of aspect ratio, displayed strong four-fold 
orientational, i.e., tetratic, correlations. Studies of the area- fraction-dependence of these 
correlations, presented in subsection 13.11 were done using the shorter rods: since the 
system size measured in terms of particle length was greater for the short particles, 
they were less sensitive to the nature of the boundary and offered better statistics. 
We observed short-ranged tetratic order, but with large correlation lengths. The high 
aspect ratio rolling pins formed nematic phases with quasi-long-range order with power- 
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Figure 5. "Phase diagram" shows region of parameter space explored and steady 
states in each region. Vertical axis represents the two different shapes of the particles 
we used. 



law decay of orientational correlations. The low aspect ratio rolling pins formed swirly 
states which we discuss in subsection 13.4.21 Common to both nematic and tetratic 
states, was the tendency to rotate globally and persistently in one direction, presumably 
in response to some asymmetry in the boundaries of the system. This distinguishes this 
phenomenon from that observed by [221 who also observe globally rotating phases, only 
the particles forming their phase are, unlike ours, intrinsically chiral. We discuss this 
in subsection 13.4.11 We also observed striped, smectic-like phases in experiments using 
rice grains. Interestingly, we have not been able to observe the stripes using bails or 
cylinders with the same aspect ratio as the rice. 

3.1. Towards the tetratic phase 

Over the range of concentration and other parameters explored, the strictly cylindrical 
particles (i.e., fiat headed, neither tapered like the "bails" that showed nematic order, 
nor spherically capped like the spherocylinders widely studied in simulations) did not 
form ordered phases. However, their isotropic phase showed tetratic, that is, four- 
fold orientational, short-range order unaccompanied by nematic correlations. The 
range of tetratic correlations increased substantially with increasing area fraction. 
The correlation function evaluated was G'4(r) =< cos[4:{9i — Oj)] >, where 9i is the 
orientation of the rods with respect to some arbitrary, pre-selected axis, i and j are 
particle indices that run from 1 to the number of particles, N. The angular brackets are 
a time average, over pairs z, j of particles separated by a distance r. For each pair, we 
calculate G4 as a function of r± and r\\ which are, respectively, the components of the 
separation vector r between the rod centres perpendicular and parallel to the cylinders' 
long axes. r±_ and r|| are distinct for each particle in a given pair, they have to be 
evaluated for each particle, rather than each particle pair, as would have been sufficient 
if G4 was evaluated as a function of the scalar distance, r. 

Both long [A = 12.6) and short {A = 5.2) cylinders showed tetratic correlations 
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Figure 6. Tetratic correlation lengtlis (scaled by rod length L) are larger for aspect 
ratio A — 5.2 than for A = 12.6. Results shown for area fractions 87% and correlations 
as function of r±_ (defined in section Rj.lll : the behaviour with respect to r[[ is similar. 

as shown in figure El The range of tetratic order is shorter (when distance is scaled by 
cyhnder length) for the long cylinders, however, even a correlation length = 1.57L = 
19.8D, as found for the long cylinders, indicates substantial tetratic ordering with clearly 
identifiable four-fold motifs as seen in FigjSJ When the experiments were performed in 
a square, rather than a circular frame, one might have expected four-fold orientational 
order to be in general enhanced, since alignment parallel to the walls is (see section 
favoured by the mechanics of collisions with the wall. We find that this is true for both 
long and short cylinders but not at all for the bails. 

To distinguish between tetratic and nematic ordering, we have plotted for the 
short cylinders A = 5.2 at area fraction 90% , G4(r) and G2{r) in Fig. d here 
G2{r) =< cos[2{9i — 9j)] > quantifies the extent of nematic ordering. G4,{r) clearly 
dominates, and nematic correlations are negligible. This confirms what is apparent 
from the real space image of the tetratic phase, figure IH We are thus seeing true 
four-fold correlations, not angular harmonics of 2-fold correlations. 

We now present the dependence of the range of tetratic correlations on area fraction 
for the short cylinders {A = 5.2). Fig. [Hlshows G4 as a function of r\\ and r±_ for various 
area fractions. The range of tetratic correlations increases monotonically with area 
fraction, as shown in Fig. ^ which displays the correlation length, ^ , extracted from 
Fig. ISl ^ obtained from r± and r\\ show the same qualitative behaviour, but are slightly 
different in magnitude. 

The correlation lengths are substantial at high area fractions, indicating a high 
tetratic susceptibility, and the likelihood that a quasi-long-ranged tetratic phase lurks 
nearby in parameter space. However, it remains to be determined which parameter - 
area fraction, confinement thickness, vibration amplitude or frequency - will most easily 
take the system across this conjectured isotropic-tetratic phase boundary. We reiterate 
that the tetratic ordering is a strong function of particle shape; by way of comparison, 
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Figure 8. Gi (defined in section IXT)) as a function of and r\\ for different area 
fractions, for cylinders with aspect ratio 5. 

the equilibrium phase diagram of [T] for discorectangles at similar aspect ratios shows 
no tendency to tetratic order. They find isotropic, nematic, and crystalline phases, 
none of which we see with rectangles. This difference could be a nonequilibrium effect 
rather than a shape effect, however as we will see in section IT2| tapered particles order 
differently. Also, Monte Carlo simulations of rectangles |2I] at much smaller aspect 
ratio, A = 2, reveal a true tetratic phase. While it is not clear at all whether this 
equilibrium phase will survive to the large A regime, we hope that our results will 
stimulate equilibrium simulations of this region of parameter space. 
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Area Fraction (%) 



Figure 9. Tetratic correlation lengtlis increase with area fraction, suggesting that 
a tetratic phase should appear if some other parameter such as shaking amplitude is 
varied. We do not know the origin of the slight difference between the correlation 
lengths obtained along r±^ and r\ [ . 



Table 2. The table below lists the correlation lengths ^ (scaled by D, the rod 
thickness) and decay exponents a for, respectively, the short-range-ordered and power- 
law-ordered nematic phases. 



Area Fraction 




a(r||) a(r^) 


50% 
60% 
70% 


69 64 


0.0113 ± 4e-4 0.0182 ± 8e-4 
0.0066 ± 3e-4 0.0140 ± 5e-4 



3.2. The nematic phase 

Why do cylinders not show nematic order even when they are extremely long? Nematic 
order is typically investigated in simulations using sphero cylinders, discorectangles or 
ellipsoids, all of which have narrower tips than waists and can thus interdigitate, 
promoting uniaxial alignment. This indicates that the rod tips needed to be tapered 
to obtain a nematic phase. We therefore used the particles labelled "bails" in figure HI 
with aspect ratio A = 10.4, since equilibrium Monte Carlo studies tell us that ^ > 7 is 
needed to produce a two-dimensional nematic. 

As is evident from the real space images (Fig. the long bails do indeed 
show nematic correlations. The tendency to nematic order is quantified in Fig. ^| 
which shows G2 as a function of r± and r\\ for 3 different area fractions 0. Here 
G2{r) =< cos[2{9i — 9j)] >. The symbols are the same as defined in section ITTl The 
system shows a high degree of order at each of these area fractions. The decay of G2 
was consistent with a power-law for all three values of (p, although the decay exponents 
(listed in Table |2)) were exceedingly small for the two higher values of 0. 

Given the small system sizes in our experiment, we expect strong finite size effects. 
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Figure 10. G2{r) (defined in section as a function of r|| and r±: The phases are 
very highly ordered and only a detailed subsystem analysis confirmed the quasi-long- 
ranged (rather than true-long-ranged) order at area fraction — 60% and 70% and 
short-ranged correlations at </> = 50%. The decay exponents and correlation lengths 
extracted from the plots are summarised in tabled 
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Figure 11. Binder reduced cumulant Ul vs area fraction suggests quasi-long-range 
nematic order at and above (p — 60%, short-ranged order below. Here L is the linear 
size of the subsystem in units of the rod diameter d 

In order to assess more precisely whether the correlations were long-ranged, power-law, 
or short-ranged with very long correlation lengths, we performed a subsystem analysis 
of the cumulants following j^Ij. The quantities evaluated were S =< ip >l, where ip 
is the order parameter {ip = e*^^^ , where j runs over the number of particles 

within the subsystem defined by L) and L is the linear size of the subsystem over which 
S was evaluated, S*^ =< ip"^ >l and S*^ =< ip'^ >l- The Binder reduced cumulant, 
Ul = 1 — g^g2^2 was then evaluated. The results are plotted as a function of area 
fraction, 0, for different subsystem sizes L in figure ^2 

As shown in 12^, for phases with LRO the curves for different system sizes L will 
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Figure 12. Binder reduced cumulant Ul linear size L of subsystem over which Ul 
was evaluated. For </) = 50 %, mild dip in curve suggests short-ranged order. Inset: 
Ul for (p ~ 60% and 70% at a higher resolution, (p is very slightly, but systematically 
higher for (p — 70%, leaving open the possibility of true long-range nematic order in 
this driven two-dimensional systems. Here L is the linear size of the subsystem in units 
of the rod diameter d 



cross at a certain value of area fraction (p = (jf . If there is QLRO, the curves of Ul 
vs for different L should collapse onto a single line over an extended range of values 
of 0. The lower limit of this value will be 0*, the isotropic- nematic transition value. 
Figure ^2 suggests that the phases formed at the higher area fractions do show QLRO 
and the isotropic- nematic transition occurs at an area fraction somewhere between 50% 
and 60%. In Fig^Jwe show the Binder cumulant as a function of subsystem sizes L, 
for different area fractions, 0. The downturn at large L for = 50% indicates short- 
range order. We show a greatly magnified version of this plot for the two higher area 
fractions as an inset. Ul for 70% is slightly but systematically larger than that for 
60%. This difference is at the edge of our resolution, but it leaves open the tempting 
possibility that these nonequilibrium nematics have true, rather than quasi- long-range 
order. Further measurements on large systems are required to settle this important 
question of principle. 

3.2.1. Varying amplitude of oscillation The presence of a third dimension into which 
the particle is lofted means that the 2-dimensional projections of the particles can 
overlap and that unlike hard particles, a temperature variable could be relevant to 
the order. As we said in subsection 12. 2| the shaking amplitude Oo could act as a kind 
of temperature. A natural dimensionless measure of how strongly the system is being 
shaken is F = {27ifYao/g, where / and g are, respectively, the oscillation frequency and 
the acceleration due to gravity. We studied the effect on the nematic order of changing 
acceleration while holding fixed the gap height, H and particle diameter, D; changing 
these would have changed the inter-penetrability of the particles in two dimensions. 
Figure [121 shows the Binder reduced cumulant of the phase as a function of 
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Figure 13. Binder reduced cumulant Ul versus dimensionless acceleration. The loss 
of long-ranged order upon increasing amplitude confirms the notion that a temperature 
variable is relevant in this system and is related to the shaking strength. 



dimensionless acceleration. The area fraction of the phase is about 63%, within the 
nematic regime identified in FiglT^ The acceleration is varied by holding the frequency, 
/, fixed and varying the amplitude, Oq. Even though << {H — D) over the entire 
range of amplitudes, the higher shaking amplitudes do affect the extent of order in the 
system. This is not entirely surprising: while Oq is small, the ratio of the characteristic 
kinetic energy l/2m(27r/ao)^ and the maximum possible potential energy mg{H — D) 
varies from 0.2 to 0.36 over the range of the plot. However, as can be seen by the small 
range of vertical axis, the degree to which the phase disorders by increasing T from 5 to 
8 is much less than that achieved by reducing the area fraction from 70% to 50% . 

3.3. Two-dimensional smectic order? 

Our observations on basmati rice gave images with strong indications of striped, i.e., 
smectic order (see Fig. In two-dimensional systems at thermal equilibrium, we know 
that thermal fluctuations render smectic order strictly short-ranged, while in three- 
dimensional systems it is quasi-long-ranged. These conclusions are based on derivations 
of the variance of layer displacement fluctuations from an elastic free-energy functional, 
in which, as a consequence of rotation-invariance, the cost of layer displacements with 
wavevector q parallel to the layers varies as g^. Since the systems we are studying 
are driven, not thermal equilibrium, these arguments do not apply a priori, and the 
possibility of quasi-long-range (or, for that matter, true long-range) smectic order in 
our vibrated granular- rod monolayers remains open. We are studying this possibility 
theoretically [28|. Experimentally, however, our systems are not large enough for the 
type of analysis (see, e.g., [201130]) required to check whether the images we are seeing 
are true smectics, weakly ordered two-dimensional crystals, or smectics disordered by 
the presence a finite but small density of free dislocations (which would be nematics on 
large length scales). 
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3.4- Dynamical nonequilibrium phenomena 

3.4-1 ■ Collective, coherent rotation All the phases showed a marked tendency to rotate 
globally. For the nematic phases, the square frame prevented them from rotating (except 
in certain cases, which we shall come to in the next paragraph) , but we could not prevent 
the tetratic phase from rotating, in part because we did not use a square frame, since we 
wanted to see how strong the tetratic correlations were in the absence of such externally 
imposed four-fold alignment. 

Since this is a nonequilibrium system, any stray chirality in the system, even 
at the boundary, should result [HI] in macroscopic circulation. The observed global 
rotation is thus a signature of the driven nature of the system. A thermal equilibrium 
system, no matter how chirally asymmetric, would not display persistent currents. Our 
observed rotation rates were of the order of a degree in 10^ cycles of oscillation. In one 
instance, when we were using the bails in the square boundary (where, under normal 
circumstances, the rotations are curbed), the phase began rotating fast, turning through 
a degree in ?5i 10^ cycles of oscillation. The source of the asymmetry appeared to be 
a slight misalignment of the particles fixed at the boundary, because we were able to 
stop the rotation by aligning them more carefully. That this minor imperfection in 
boundary alignment could have such a strong effect was a surprise, and means the 
rotation phenomenon merits more serious investigation. 

3.4-2. Swirls Since basmati rice (aspect ratio A = 4.2) showed smectic tendencies, we 
studied tapered copper rods with similar aspect ratio (aspect ratio A = 4.9) but half 
the length of the rice grains, in the hope of getting a larger number of smectic layers 
in the sample cell. Unexpectedly, we found no sign of smectic layering. Instead, a very 
dynamic pattern of swirls appeared, reminiscent of bacterial swimming patterns |32j . 
The swirly state looks like a defect- ridden nematic, which is curious because equilibrium 
Monte Carlo studies on rods in two dimensions suggest that nematics should not form 
for aspect ratios less than 7. These swirls, which are strength ±1/2 disclinations, differ 
from those seen in P^: the latter are strength 1, they have a core that looks like a 
crystalline packing of standing rods, and each defect is chiral because the rods have 
a systematic tilt in the direction of their tangential velocity orientation of the rods in 
the region outside the core. In our case, the circulation appears to be a result of the 
interaction between defects. It remains to be seen whether we are seeing coarsening 
towards a nematic, the disordering of a nematic (i.e., coarsening towards the isotropic 
phase), or simply a disordered, correlated, defect-turbulent steady state. 

4. Conclusion and Discussion 

This paper has presented an exploration of the collective properties of vibrated rods 
in two dimensions, focusing on orientationally ordered phases. Tapered particles lead 
to nematic order at large area fractions, in keeping with the thermal equilibrium 
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Monte Carlo studies of P, while strictly cylindrical particles, with flat tips, display 
strong tetratic correlations over a very broad range of aspect ratios. It remains to 
be established whether the strong sensitivity to particle shape is a nonequilibrium 
phenomenon or whether it can be reproduced in equilibrium simulations. We also 
observed some clearly nonequilibrium phenomena: persistent, systematic rotation of 
the entire phase, presumably as a result of stray chirality in the system, and also a rich 
pattern of circulating disclinations of strength ±1/2. Detailed quantitative analyses of 
the dynamics of the swirls, of tagged-particle motion and time-dependent correlations of 
collective modes and other nonequilibrium effects will appear in a separate paper 
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